Diffusion of a liquid nanoparticle on a disordered substrate. 
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We perform molecular dynamic simulations of liquid nanoparticles deposited on a disordered 
substrate. The motion of the nanoparticle is characterised by a 'stick and roll' diffusive process. 
Long simulation times (~ /xs), analysis of mean square displacements and stacking time distribution 
functions demonstrate that the nanoparticle undergoes a normal diffusion in spite of long sticking 
times. We propose a phenomenological model for the size and temperature dependence of the 
diffusion coefficient in which the activation energy scales as N^'^. 
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I. INTRODUCTION 

The spatial organisation of nanoparticles deposited 
on a substrate strongly depends on the nature and the 
magnitude of their mobilty during the growth process. 
For systems known as Volmer- Weber films [l| for which 
nanoparticles are isolated from each other, it is of im- 
portance to understand their diffusion on the substrate. 
This diffusion is directly related to the interaction be- 
tween the substrate and the nanoparticle so that depend- 
ing on the system studied one can found slow or fast dif- 
fusing nanoparticles la, |3, II] . This has been illustrated 
by Deltour et al. ^5' who demonstrated that a small 
change in the mismatch between the lattice parameters 
of a Lennard- Jones cluster and the substrate induces an 
important variation of the diffusion coefficient. More re- 
cently, a surprising increasing mobility with increasing 
size of the nanoparticle has also been reported [g 

From a more fundamental point of view it is impor- 
tant to determine the nature of the nanoparticle diffusion. 
The simplest method is to look at the time dependence 
of the mean square displacement of the center of mass : 
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where Ycm{t) is the position of the center of mass at time 
t. The exponent 7 is equal to one for a normal diffu- 
sion while 7 > 1 and 7 < 1 correspond respectively 
to superdiffusive and subdiffusive processes. Luedtke 
and Landman J3, Q have performed simulations of gold 
nanocrystals adsorbed on graphite and exhibiting Levy 
Flights leading to mild superdiffusion (7 ~ 1.1). More re- 
cently and for a similar system Maruyama and Murakami 
[3 found a crossover between superdiffusion and normal 
diffusion. These different diffusion regimes are influenced 
by the way the nanoparticle moves on the surface. Slid- 
ing of the whole nanoparticle seems the most probable 
scenario but a rolling mechanism has also been observed 
for small clusters 
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FIG. 1: Snapshot of a TV = 555 gold liquid nanoparticle de- 
posited on the disordered substrate. We measure a wetting 
angle 6 ~ 140° 
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To our knowledge and in contrast to the case of solid 
nanoparticles, a single study [ij has been devoted to 
the diffusion of liquid nanoparticles. The authors have 
reported a discontinuity in the diffusion coefficient ob- 
tained for solid and liquid nanopaticles near the melt- 
ing point. In this paper we perform molecular dynamics 
simulations of gold liquid nanoparticles deposited on a 
weakly disordered substrate. The first section of this pa- 
per presents the simulation details and in the second one 
we characterize the motion of the nanoparticle on the 
substrate. We show that it is made of sticking and mov- 
ing events. The existence of a contact angle hysteresis 
is at the origin of the oscillations around an equilibrium 
position. A visual inpection and the measure of the cor- 
relation between the velocity and the angular momentum 
evidence the tendency of the nanoparticle to experience 
a 'rolling like' motion to reach another equilibrium posi- 
tion. In the third section we examine the nature of the 
diffusion. Long simulation times are necessary to com- 
pute reliable mean squared displacements and sticking 
time distribution functions. We show that the diffusion 
remains normal even for the largest nanoparticles having 
large sticking characteristic times. A phenomenological 
model is finally presented to account for the size and 
temperature dependence of the diffusion coefficient. The 
main hypothesis of a size dependent activation energy is 
necessary to recover the simulation results. 
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FIG. 2: Time evolution of the x coordinate of the center of 
mass for a particle of 555 atoms. The sticking and moving 
events are clearly identified. Inset : during a sticking event 
the CM oscillates around an equilibrium position with a char- 
acteristic vibrational period r„ . 



II. SIMULATION PROCEDURE : 

The nanoparticles considered in this study have sizes 
ranging between N = 87 and N = 555 atoms. We use the 
Ercolessi glue potential ^^ to describe the interactions 
between Au atoms with a cut-off distance of 3.9A and a 
time step St = 2.5 lO^^ps. In this formahsm the potential 
energy of the N atoms is given by : 
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The first term is a classical pair interaction. In the second 
term, rii is the coordination of atom i : 



N 



= E^(^*j) 



(3) 



where p{rij) is a function of the interatomic distance rij. 
The energy function U is the glue term associating an 
extra potential energy which is a function of the atom 
coordination. 

The substrate is a slab of dimension 170 * 170 * lOA'^ 
made of 15562 atoms. The substrate atomic positions 
come from the simulation of a liquid rapidly quenched 
below its melting temperature. We obtain a substrate 
with a glass structure and periodically boundary con- 
ditions applied in the x and y directions. The substrate 
atoms are chosen to be frozen in order to increase the sim- 
ulation time. We therefore expect a quantitative change 
of the diffusion coefficient as compared to a more real- 
istic simulation in which the substrate atoms can move 
[l3| . Nevertheless this choice permits to reach simulation 
times of roughly 1/is that, as will be discussed below, are 
necessary to clearly identify the diffusion regime. The 



interaction between substrate and Au atoms are chosen 
to be of the Lennard jones type : 
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where r is the distance between the two atoms. We 
choose a — 2.7 A, e — 255 fcf, and a cut-off distance of 
2.5 cr. As a consequence of this weak e value, the liquid 
Au nanoparticles are in a weakly wetting situation with 
a wetting angle roughly estimated to 9e — 140°. A snap- 
shot of the simulated system is given in Figure 1 . where 
we can see a A^ = 555 liquid nanoparticle deposited on 
the disordered substrate. The equations of motion are in- 
tegrated using the Verlet algorithm and the temperature 
of the nanoparticle is fixed using the classical velocities 
rescaling procedure [l3|. We choose T = 800 JiT below 
the bulk melting point (T,„ ~ 1380 A') but still above 
the solidification temperature. The liquid state of the 
nanoparticles is evidenced by the atomic mean squared 
displacements calculated in the frame of the center of 
mass. The values obtained for the atomic diffusion are 
in good agreement with previous ones |l2|. For exam- 
ple, for a particle of 249 atoms the diffusion coefficient 
is estimated to be 0.09 A^ps~^. In a first run of roughly 
10^ MD steps the liquid nanoparticle is deposited on the 
substrate and reach its equilibrium shape. A second long 
run of 10^ MD steps is then performed. It permits to 
record the position of the center of mass and to compute 
different quantities like the probability density function 
of the sticking events. We checked that the same results 
are obtained if the second run is performed under micro- 
canonical conditions. The velocity rescaling procedure 
has therefore no infiuence on the nanoparticle motion. 



III. ANALYSIS OF THE LIQUID DROP 
MOTION 

The trajectories of the center of mass (CM) show that 
the diffusion process is composed of two distincts events 
: sticking events in which the CM oscillates around an 
equilibrium position and moving events in which the CM 
reaches another equilibrium position. This is illustrated 
in fig. 2. where we plot the time evolution of the CM x 
coordinates Xcm{t) recorded during 20ns for a particle of 
555 atoms. For this large particle, the sticking and mov- 
ing events are clearly identified. In inset we represent an 
hundredfold time expansion of Xcm{t) for the same sys- 
tem. During a sticking event the CM oscillates around 
an equilibrium position with a characteristic vibrational 
period Ty. A Fourier transform oi Xcm(t) and ycm{t) per- 
mits to evaluate t„. These results will be presented else- 
where but as a main result the vibrational period t^ is 
found to scale with the square root of the particle size. 
Indeed because of the hysteresis contact angle [lj| the 
liquid drop is at rest on the substrate and is exited by 
the thermal noise. The measured t^ corresponds to the 
eigen vibrational period due to the elastic deformation of 



the drop : 
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(5) 



where p is the density and 7/^, the hquid-vapor surface 
tension. This period is similar to the one found for 
the shape oscillations of a free liquid drop [13 • When 
the particle leaves its equilibrium position, we observe 
a 'rolling like' move (note that the liquid drop is also 
changing its shape during the move) of the liquid drop 
on the substrate. In fig. 3 we represent three snapshots 
of a moving particle at times separeted by 50 ps. The 
atoms initially belonging to the lower half of the particle 
are grey colored to evidence the 'rolling like' motion. To 
quantify this visual inspection we calculate the correla- 
tion between the components of the CM velocity {Vx and 
Vy) and the orbital moments {Jx and Jy). As usually the 
correlation coefficient C^ between the components i and 
j are defined as : 
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FIG. 4: Mean square displacement -Rcm(i) for particles with 
sizes A*' = 87, 249 and 555. The straigth lines correspond to 
the asymptotic linear regimes : Remit) c< 4Dt. 
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We report on Table I the different C^ coefficients com- 
puted for a, N = 321 particle during a run of 150ns. 
The correlation is not significative between Vx and Jx 
and between Vy and Jy. Conversely and as expected 
for a 'rolling like' move, a correlation is found for the 



cross components of the velocity and the orbital moment 
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FIG. 3: Snapshots of a moving nanoparticle taken at times 
separeted by 50ps. The atoms initially belonging to the lower 
half of the particle are grey colored to evidence the 'rolling 
like' move. 



Note finally that these correlations coefficients are 
surestimated because they also characterize the vibra- 
tions around the equilibrium position. Therefore and 
since the radius of our particles should be of the same 
order than the slipping length [ifiL \u\ , we cannot rule 
out possible sliding moves of the nanoparticle. 



IV. SIZE AND TEMPERATURE DEPENDENCE 
OF THE DIFFUSION COEFFICIENT 

We now present the results obtained for the mean 
square displacement i?^„(i) of particles with sizes rang- 
ing from TV = 87 to TV = 555 atoms. We plot in figure 
4. the statistically rehable values of i?^„(i) ( i < 20 ns) 
computed for a run length of roughly 1/is. For the small- 
est particle (N = 87) the regime of normal diffusion is 
clearly indentifiable with B?{t) oc t. As the particle size 
increases the time to reach a linear regime increases. For 
the largest particle considered in this study [N = 555) 
the linear regime is not reached and the conclusion from 
the B?{t) curve is not clear : one could think at a subd- 
iffusive regime with a fitted exponent 7 ~ 0.72 or alter- 
natively that the simulation time is not long enough to 
evidence a normal diffusion regime. 

In order to reach a conclusion about the nature of the 
diffusion for the largest particles, we look at the probabil- 
ity density function (PDF) pr{t) of the sticking events. 
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FIG. 5: Rank ordering of the sticking times for a particle 
of 249 atoms {Ra = 25 A). The slope of the straight line 
corresponds to the mean sticking time of the exponential dis- 
tribution. Inset : the diffusion constants D{Rs) tend to a 
constant for sufficiently large Rs values. 



Indeed we know [1^ that the diffusion nature depends on 
the form of pr{t) : for a PDF that has a mean value (like 
an exponential one) the diffusion is known to be normal 
while for a power law PDF ( oc r^(^+^' with < /i < 2) ) 
we expect a subdiffusive behavior. The time r^ of a stick- 
ing event is defined as the necessary time for the CM to 
cross a circle of radius Rg centered on its initial position. 
We compute the PDF for different particle sizes and dif- 
ferent values of Rg ■ Since the statistic on sticking times 
is rather poor especially for the largest systems, the na- 
ture of the PDF is determined using the rank ordering 
method J19l |. For an exponential distribution we expect 
that the the n — th largest value of t^ , t" scale as : 



-fslog(n/N) 



(7) 



where N is the total number of observed sticking events 
and fs the mean sticking time. As can be seen in figure 
5. for a particle of 249 atoms and Rg = 25A the func- 
tion r" is clearly logarithmic in n. Whatever the particle 
size considered, for sufficiently large Rg values the asso- 
ciated PDF is exponential and we can measure a mean 
sticking time f^ . This analysis demonstrates that the dif- 
fusion remains normal even for the largest particles. The 
apparent subdiffusive behavior observed in the i?^^(i) 
curves is due to the large value of the mean sticking time 
as compared to the simulation time. We also compute 
the mean square displacement f^ between two sticking 
events and finally estimate the diffusion coefficient to be 



D{Rs) 



' I^Tg. We can see on the inset of figure 5. 



that D{Rs) tends to a constant for sufficiently large Rg 
values. We establish here a first method to measure the 
diffusion coefficient. 

Another possible way to extract D{N) from the sim- 
ulation data is to look at the asymptotyc behavior of 



R^^{t). We have demonstrate above that the diffusion 
remains normal whatever the particle size so that the 
quantity 4>(t) = i?^„j(i)/4i should be of the form : 



$(i) = D + ar 



(8) 



for large t values with an exponent 7 < to ensure the 
asymptotic normal diffusion regime. Best fits of our data 
to this form of <l>(i) give values of 7 close to —0.5 with 
a small dispersion around this mean value (cr^ ~ 0.05). 
This exponent is the same than the one predicted for 
the model of continuous time random walk ll8| with 
a sticking time distributed according to a power law 
p{t) oc T^^^+t^^ with /i < 3 (normal diffusion regime). We 
therefore extract the diffusion coefficients D{N) through 
best fits of simulated ^(i) values to equation 8 (the value 
of 7 is fixed to 0.5 for the fitting procedure). We repre- 
sent in figure 6. the simulated and fitted ^{t) values for 
systems with N = 177, 249 and 429 atoms. The arrows 
indicate the values of the diffusion coefficient extracted 
from this asymptotic analysis. 

We represent in figure 7. the diffusion coefficients ob- 
tained by the analysis of the sticking times PDF and by 
the second method described just above as a function 
of the system size. A first fit of these data to a power 
law D{N) ex A^T leads to a value 7 = -1.3 ± 0.1. The 
agreement is satisfactory for the smallest particles but we 
observe a deviation for the largest N values. The power 
law apparently tends to overestimate the value of D. For 
a LJ cristalline nanoparticle on a solid substrate Q the 
exponent is 7 = 2/3 for a large mismatch parameter. 
This exponent has been recently confirmed for gold sup- 
ported nanoparticles |23|. In the case of quasi epitaxy 
the exponent increases apparently up to ~ 1.4 Q. This 
latter value is close to the one found in the present study 
for a liquid nanoparticle. In both cases, a critical applied 
force exists, below which the system cannot move. In 



0,06 



^" 0,04 



0,02 




e 



t(ns) 

FIG. 6: Assymptotic analysis of >l>(f) = Rl,,-,{i)/4i for 
nanoparticles with N = 177, 249and429. Full lines are best 
fits to the simulated data and the arrows indicate the extrap- 
olated values of D. 



our case this force is related to the wetting contact angle 
hysteresis. To our knowledge no theoretical prediction 
of this large exponent value exists. Note finally that the 
diffusion coefficients calculated with the fitted power law 
for large particle sizes ( for example Z? ~ 4 nm?s~^ for 
a particle of radius 2± IQQnm ) are definitely too large 
to be consistent with experimental observations J2ll |2^ 
of completely sticked liquid nanoparticles. This suggests 
that the power law is no longer vadid for large particles. 
We therefore propose a second phenomenological scal- 
ing law of the form : D{N) oc exp{—aN^/^/kT). In this 
expression the activation energy of the diffusive process is 
a function of N . To escape from its equilibrium position 
the liquid drop is elastically deformed with advancing 
and receding wetting angles differents of the equilibrium 
one. The exponent f/3 in the proposed scaling law re- 
lates the elastic deformation caused by the wetting angle 
hysteresis and the parameter a is therefore proportional 
to 7ii,. We can see on fig. 7. that this scaling law reason- 
ably fit our data on the whole range of particle sizes with 
a free parameter a = {).QbleVatoms~^/^. As done just 
above we can estimate the value of D for large particles. 
This time the result found is in agreement with the ex- 
perimental observation |2lL l2a | . Note that we have made 
the hypothesis that the main N dependence is in the 
exponential. Nevertheless a size and temperature depen- 
dent prefactor should also exists so that a more general 
expression should be of the form : 



D{N,T) ex e-^^-A^/(iV)e-'^^=W/^^ 



(9) 



The first exponential term contains the activation en- 
ergy of individual atoms. In the present case of a liq- 
uid nanoparticlc this activation energy is related to the 
atomic diffusion within the liquid. The two other terms 



have been discussed above, they are reflecting the collec- 
tive diffusion of the nanoparticlc. The form of the func- 
tion f{N) depends on the nature of the dissipation and 
could be evaluated in the same way than for a millimetric 
drop PI. 

To confirm this exponential relation between the dif- 
fusion coefhcient and the particle size we finally look at 
the temperature dependence of D for a particle of 177 
atoms. The diffusion coefficients are computed for tem- 
peratures between 700 and I200iir. Data are well fitted 
by an arhenius law D ex exp{—AE/kT) and give an ac- 
tivation energy AE = 0.504ey (Fig. 8.). This value is 
greater than the one obtained from the previous analysis 
of the D{N) data. Using the fitted a value we indeed 
found AEc{N = 177) = 0.286el^. Equation 9. predicts 
that the overall activation energy AE should be the sum 
of AEi and AEc- We then perform simulation to mea- 
sure the atomic diffusion Di within the liquid. The values 
Di (T) are extracted from mean square displacement cal- 
culated in the frame of the CM for different temperatures. 
An arhenius law fits the data and gives the value of the 
individual activation energy : AEi = 0.196eT^. Finally 
we found AEc{N = 177) + AEi = 0.482ey that compares 
rather well with AE = 0.504ey and seems to confirm the 
phenomenological form proposed for D{N, T) in Equ. 9. 



V. CONCLUSIONS AND PERSPECTIVES 

To summarize, we have performed molecular dynamic 
simulations of a liquid nanoparticlc diffusing on a weakly 
disordered substrate. For the weakly wetting condition 
considered here, the nanoparticlc experiences a 'rolling 
like' motion from an equilibrium position to another one. 
The sticking time distribution function is shown to be 
exponential even for the largest nanoparticles. As a con- 
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FIG. 7: Diffusion coefficients as a function of the system 
size. We represent the value obtained by the two methods. 
The two straight lines correspond to the best fits to a power 
law and to a law in which the activation energy is a function 
of the system size. 
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FIG. 8; Diffusion coefficients as a function of temperature 
for a, N = 177 particle. The best fit to an arhenius law gives 
an activation energy of 0.504eV. 



sequence the diffusion is normal with a mean square dis- 
placement asymptotically linear in t. The diffusion coeffi- 
cients are reported as a function of size and temperature. 
A phenomenologicalexpression for D{N,T) is proposed. 
The main hypothesis of an activation energy scaling with 
N^/^ is necessary to fit the simulation results. 

In spite of the difficulty of a direct observation of the 
nanoparticle diffusion, we hope this paper will motivate 
further experimental studies. Finally, a work is currently 



in progress to fully characterize the vibrational modes of 
the sticked liquid nanoparticle. 
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